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Abstract. In this work we apply field regularization techniques to formulate a number of new phenom- 
ena related to momentum induced by electromagnetic zero-point fluctuations. We discuss the zero-point 
momentum associated with magneto-electric media, with moving media, and with magneto-chiral media. 
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1 Introduction 

It is well-known and widely accepted that zero-point fluc- 
tuations affect the physics on both microscopic and macro- 
scopic scales, upon creating forces between materials. The 
Casimir effect [1|- the attractive force between two macro- 
scopic metallic plates - and physically equivalent to the 
Lifshitz effect when it comes to dielectric media - is un- 
doubtedly the most famous effect. Also the 1/r 6 Van der 
Waals force and its retarded 1/r 7 equivalent, the Casimir- 
Polder force [2], between microscopic polarizable atoms 
can be understood as direct manifestations of zero-point 
energy [3]. 

Electromagnetic zero-point energy has been the sub- 
ject of many fundamental work, sometimes heavily de- 
bated in literature, since it often lacks experimental ver- 
ification. The most fundamental aspect of vacuum fluc- 
tuations is their Lorentz invariance. It can be shown that 
only isotropic radiation with power spectrum w 3 , obviously 
the one associated with modes whose density is propor- 
tional to w 2 /cq and whose zero-point energy is \hu> [4] 
is Lorentz-invariant. The Unruh effect [6] is another fun- 
damental result and concerns observers with constant ac- 
celeration a. They see the zero-point energy emerge as a 
Planck law with temperature kT = ha/2-KCQ 011] ■ This 
effect has never been observed. On the basis of the rel- 
ativistic equivalence principle, gravity should create the 
same effect, known as Hawking radiation[7]. 

The most controversial aspect of vacuum energy is the 
UV catastrophe [5]. The vacuum energy density is arbi- 
trarily large at large frequencies. No rigorous mathemat- 
ical tool seems to exist so far to deal with this problem. 
Fortunately, the divergence does not affect the Casimir 
force, since it formally drops out, as can be seen for in- 
stance using the Euler-Maclaurin summation formula 0]. 
Yet, the energy itself may still be an observable. In 1993, 
Schwinger [5] attempted to explain the electromagnetic 



energy observed from strongly oscillating water bubbles 
- an acousto-optical effect called sonoluminescence - in 
terms of the zero-point energy released by the contract- 
ing bubble. Upon disregarding vacuum fluctuations with 
frequencies beyond the UV, Schwinger concluded that the 
excess Casimir energy of a sphere with radius a and di- 
electric constant e be equal to E c ~ ha?Lo^/c^{l — l/\/s) 
which has the right scaling - it decreases with decreasing 
volume - and about the right order of magnitude if the cut- 
off is properly chosen. Another longstanding issue, first 
raised by Dirac in 1934 [10], is the observation that zero- 
point energy should be gravitationally active and should 
thus appear as a contribution to the cosmological con- 
stant in the Einstein equations. If this is true, the diver- 
gence of zero-point energy comes in rudely. The Planck 
length ^/TiG/cl = KT 33 m is the only available cut-off for 
vacuum modes hand but this would still lead to impossi- 
ble cosmological scenarios. Dirac concluded that the large 
zero-point energy is absent for a still mysterious reason. 

So good motivations exist to search for regularization 
techniques that deal with the divergence in a different 
way. Milton et al [11 proposed a more sophisticated reg- 
ularization scheme, where only the finite part of the zero- 
point energy is considered, due to the finite geometry [13] . 
For a spherical bubble they find E c ~ +h(e — l) 2 /aco, 
which has clearly the wrong scaling to explain sonolumi- 
nescence, and which is orders of magnitude smaller than 
the estimate by Schwinger. The same regularizations have 
been tested for a variety of space-time structures to solve 
the cosmological constant problem (see [TJ] for a review) . 
Here the sign is also an issue since the cosmological con- 
stant is believed to be positive, leading to a negative pres- 
sure and an expansion that accelerates. The sign of the 
Casimir force itself had already been an issue in the six- 
ties. In 1956 Casimir himself proposed [T5j that the vac- 
uum force exerted on a spherical metallic shell with sur- 
face A = Aira 2 might have the similar attractive form 
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F = —ahcoA/^ira 4 — —ahco/a 2 as was found for the 
plates, though with a different unknown constant a > 0. 
He speculated that this force might stabilize the Coulomb 
repulsion F = +e 2 /a 2 of the electron. This would provide 
a first calculation for the fine structure constant, since 
stability would impose e 2 /hcQ = a. Unfortunately, the 
regularized Casimir force on a metallic surface was shown 
by Boyer to be repulsive [5] . 

Only a few years ago, a new controversial effect trig- 
gered by zero-point fluctuations was put forward, this time 
addressing their momentum in magneto-electric (ME) me- 
dia. In any medium the dielectric constant can be af- 
fected by external electric and magnetic fields according 
to As — • (Eo x Bo), with k the unit wave vector. This 
leads to different optical properties for photons propagat- 
ing along or opposite to the vector E x B = So though, 
unlike the Faraday effect, independent on circular polar- 
isation. Rcf. [16 considered the radiative momentum of 
photons in ME media. His predictions are made in a con- 
text that is already controversial in itself, since the mo- 
mentum of photons in matter is still heavily debated [17] . 
The momentum density of zero-point fluctuations in a ME 
medium was found to be, 

Like Schwinger in his attempt to explain sonolumines- 
cence, Feigel regularized by adopting an UV cut-off for the 
zero-point spectrum, arguing that at very high frequen- 
cies the ME optical response should vanish. His choice of 
a lower cut-off wavelength of 0.1 nm is based on the fact 
that optical ME has been observed in the X-ray. He as- 
sumed a typical ME effect x^o ~ 10 -11 , and with mass 
densities typically equal to 1 g/cm 3 this would lead to 
typical speeds of v = 30 nm/s, likely to be too small to be 
measurable. However, our literature study revealed that 
ME materials exist such as FeGa03 for which xSq ~ 10 -4 
is observed down to wavelengths of order 2 A in the X- 
ray [18] . The prediction in Eq. ([1]) would lead to much 
larger speeds, up to centimeters per second, that should 
be observable in experiments. 

The calculation of zero-point momentum in ME media 
puts forward a revolutionary prediction (by APS Focus 
|19| referred to as "momentum from nothing" ), with a 
clear order of magnitude estimate. This work is a new oc- 
casion to question cut-off procedures for zero-point modes. 
They break the Lorentz-invariance of the quantum vac- 
uum provocatively and indeed the end-result (1) is so 
much Lorentz- variant, that it is not even likely to be re- 
pairable. But most of all, like in the Schwinger theory 
of sonoluminescence and in the cosmological constant de- 
bate, the cut-off procedures give "inelegant" and "unrea- 
sonable" results. As for zero-point momentum, the real 
QED vacuum is known to have a frequency-independent 
ME response x ~ ^ e / m t c o 1^3' so that Formula 1 pre- 
dicts a finite zero-point momentum density of "empty" 
vacuum up to 10 50 times larger than the momentum den- 
sity E x B /47rco associated with the applied fields. In 
matter, the cut-off procedure is often justified as a crude 



way of dealing with dispersion, but the results above sug- 
gest that this may not be the whole story, and that in 
reality the UV catastrophe is nonexistent for a yet un- 
known reason. 

In this work we investigate how the zero-point momen- 
tum emerges if one applies the field regularization tech- 
niques that have been proposed in literature. This tech- 
nique would eliminate the Schwinger theory as an explana- 
tion for sonoluminescence rTTllT2] . It is not our intention of 
this work to advocate regularization techniques. We wish 
here to come to a quantitative prediction by assuming the 
validity of these techniques. They are well defined math- 
ematically and straightforward to implement numerically, 
even in symbolic software. However, to our knowledge no- 
body has ever been able to assign the removed, diverging 
terms to the values of observable constants, as it should 
be in a good renormalizablc theory. Also experimental 
tests are rare. Brevik etal [12] and Barton [13] regularize 
the zero-point energy of a dielectric sphere with volume 
V = 47ra 3 /3 and dielectric constant e and show that this 
method is equivalent to a dimensional regularization of 
the Van der Waals energy between the atomes constitut- 
ing the sphere, 



Jv Jv { 4vr|r-rrJ + 1536™ (2j 



In section 2 we show that this regularization appears again 
in the expression for zero-point momentum. In a previ- 
ous Letter we have already applied the field regulariza- 
tion methods proposed by Kong and Ravndal [3T] for the 
ME zero-point momentum in the Casimir geometry, and 
concluded that the effect survives the regularization, but 
that its value is reduced by some 20 orders of magnitude. 
In the present work we address a genuine finite object: a 
ME sphere. This makes the regularized expression above 
subject to experimental tests, since the momentum of a 
finite object is a measurable quantity, much more than 
energy. This would be an indirect test for the regulariza- 
tion of zero-point motion in general, and such knowledge 
could be of vital importance to proceed for instance in 
the cosmological constant debate. We will also consider 
two other situations where zero-point momentum might 
show up: a moving sphere and a magneto-chiral sphere. 
Both cases reveal a surprise. For the moving sphere zero- 
point fluctuations seem to achieve a (regularized) momen- 
tum proportional to the velocity of the sphere. This would 
thus contribute to the mass of the sphere! We will show 
that the problem of a moving dipole is actually not UV 
divergent and that a precise prediction is obtained for the 
contribution of zero-point modes to the mass of a polar- 
izable atom. Finally, for a magneto-chiral sphere, we will 
present a microscopic argument why the contribution of 
zero-point motion to momentum should vanish. We hope 
that this gives deeper insight into the microscopic nature 
of chirality. 
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2 Bi-anisotropic sphere 

Starting point of our theoretical study is the set of macro- 
scopic Maxwell equations - expressed in Gaussian units - 
applied to bi-anisotropic matter [22]. Such media are de- 
scribed by a general linear "constitutive" relation between 
the macroscopic electromagnetic fields D, H, and the mi- 
croscopic fields E, B, 



with the symmetric vacuum stress tensor T> = (EiE* + 
BiB*)/4:Tr - S i:j S with £ = (E • E* +B • B*)/8tt the elec- 
tromagnetic energy density. Upon integrating Eq. §§§ far 
beyond the physical size of the object we get, 
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ExB' + pv = lim <pdS-T (10) 



D = e-E + x-B 



H 



-x 



E 



fj, 



B 



The constitutive tensors e and \x are assumed real-valued 
symmetric, the constitutive, bi-anisotropic tensor x is as ~ 
sumed real- valued. In this first work we wish to exclude 
the presence of optical dispersion and absorption. In inho- 
mogeneous media all tensors depend on the position vec- 
tor r. Time-dependence can be allowed as well provided 
the variation is much slower than the typical cycle oscil- 
lation of the electromagnetic fields, so that we can still 
work at constant frequency. The best-known case of op- 
tical bi-anisotropy is undoubtedly rotatory power, which 
can be described by the symmetric tensor Xij = 9&ij > with 
g a pseudo scalar, induced by some microscopic chirality. 
ME media can be modelled by the anti-symmetric choice 
Xij — x(E®Bj - B®E°), with x a scalar. These relations 
are combined with two Maxwell's equations applied to 
harmonic fields exp(— iujt), 



— iu'B = +co0p • E 

— iuT) — — cq4> p ■ H — 47rJ 9 

ip ■ B = 

ip ■ D = 4irp q 

This leads to the following wave equation, 



/\ . W . T i \ ■ ^ l \ 

e{r) -i—(j)p-X 

Co c 

-cj> p ■ Mr)" 1 ■ <j> P ] ■ E = 



(3) 
(4) 
(5) 
(0) 



„2 ^1 



(7) 



in terms of the hcrmitian tensor operator </> nm , p = ie n miPi- 
From this equation we can identify the interaction between 
matter and radiation, 



V(r, p) = "~~2~ [1 e(r)] + i^-<j) p ■ x T (r) - J^xW ■ < < P 



— <P P -x (r) - 1— ; 

Co c 
<j> p ■ [1 - Mr)- 1 ] • P (8) 



Upon combining the macroscopic Maxwell-equations, the 
constitutive equations and the Lorentz- force E + v/co x B, 
we can arrive at the following momentum conservation 
law, 



/ 1 



V 4-7TCO 



E x B + pv) = V • (-pw + T ) 



(9) 



To work out the flow of momentum to infinity, ex- 
pressed by the surface integral on the right, we consider a 
radiation field I(uj, k) incident on a conservative, confined 
bi-anisotropic obstacle. In the far field the electromagnetic 
polarization is orthogonal to the direction of propagation 
dS. As a result the terms EiE* and BiB* of T can be 
seen not to contribute to momentum. After some algebra 
we find that, 



lim * dS ■ T 



duj / dk in / dkout ^(w,k in ) 



dg(k in ,k out ) 



^in kout 



(11) 



Simultaneous P-T symmetry guarantees that dcr(k^ n , k ou ^, x) 
d(r(k ou t, k^ n , x) [H] • Hence, if the radiation field I(ui, k) 
is isotropic, and this is true when the obstacle is subject to 
zero-point radiation, the momentum flow to infinity van- 
ishes. In particular, no inelastic effects occur due to recoil 
effects. As a result, the total momentum 



mv 



1 



Attq 



-E x B* 



o 



is a conserved quantity. We will refer to the first term as 
the kinematic momentum and to the second term as the 
radiative momentum. The above formula for total momen- 
tum agrees with the more sophisticated theory by Nelson 
25J and we refer to this work for a more detailed discus- 
sion in relation to the Abraham-Minkowski controversy 
[TT1I21)] about which term is the real "radiation momen- 
tum" and which part constitutes the genuine momentum 
of "matter" . In particular, the contribution of radiation to 
momentum found here is equal to neither the "Abraham 



value" -J-E x H* 

47TCQ 



nor the "Minkowski value" — !— D x 

47TC0 



B*. We emphasize that for an isotropic monochromatic 
wave field scattered from a finite object, the space inte- 
gral is perfectly finite since the integrand is confined in 
and around the object. The problems will appear when 
integrating over a power spectrum that diverges itself as 
cr 3 . 

Another pertinent remark is that the momentum con- 
servation expressed by Eq. (TiT)]) continues to be valid if 
the constitutive tensors, the tensor x{ r ) m particular, are 
time-dependent. Constitutive equations with time-dependent 
coefficients can be justified when the variation of the co- 
efficients is slow compared to the variation of the fields 
themselves. This becomes a delicate issue for vacuum fluc- 
tuations that comprise all frequencies and that can thus be 
arbitrarily slow. The regularization techniques show that 
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for an object of size a, the typical frequency that con- 
tributes to the momentum equals Co/a. In cut-off proce- 
dures even higher frequencies dominate. This leaves enough 
room to turn on the external fields adiabatically. Since a 
perfect symmetry exists between the wave vectors p and 
— p if x = 0, we anticipate E x B to vanish before turn- 
ing on the external fields. The conservation law (|10| thus 
leads us to the conclusion that after having turned on the 
fields, the object achieves a velocity given by, 

rav = — [ d 3 r (0|E x BIO) (12) 

47TC J 

where (0| ■ ■ ■ |0) stands for vacuum expectation. For the 
product of two electric fields this expectation can be ob- 
tained from the fluctuation-dissipation theorem, which at 
zero temperature takes the form, 

(0|^( W ,r)£>',r')|0) 

= - 4ft^Im G i:j (uj, r, r') x 2ttS(uj - uj') (13) 
c o 

with G the (classical) Green's tensor associated with the 
wave equation l[7]). It can be straightforwardly verified that 
the momentum of zero-point fluctuations is expressed as, 

P racM = -r^— [d 3 v (0|E x B|0}< = — ^ X 

47TC0 J TTCg 

f°° f d 3 k 
ImJ da; J uj (hGjj (w, k, k) - kjG^ {uj, k, k)) 

(14) 

For a genuine empty vacuum this reduces to the familiar 
expression 2 x J e? 3 k/(27r) 3 ^Kk, which is zero in view of 
the perfect symmetry in k. 

In the following we shall consider a sphere with a di- 
electric constant e slightly different from one and a weak 
bi-anisotropic tensor x, both confined and constant in the 
sphere. We shall expand the Green's function in the po- 
tential interaction l|8|). Only contributions linear in \ can 
survive the symmetry between k and — k. We leave tech- 
nical details to the Appendices. The first order Born ap- 
proximation to G involves one scattering from the sphere 
expressed by G«(k,k) = G («,k) • V(w,k) • G («,k) 
in terms of the free-space propagator Go(w,k). The fre- 
quency integral can easily be performed and we find, 

mvi ~ tijkXjk ° 3 J d 3 k hk (15) 

Dimensional regularization puts the integral to zero, and 
no zero-point momentum - proportional to the volume of 
the object - is found in this order. Note that this contribu- 
tion is considered by Feigel, and handled using a cut-off. 

We proceed with the second order Born approxima- 
tion, and collect the terms proportional to (e — l)x- This 



involves one normal and one bi-anisotropic scattering. We 
write P r J ad = — Ijij) with 

1 j f°° , f d 3 k d 3 k' 
[G (k) • V(w, k, k') • Go(k') • V(w, k', k) • G (k)] j7 

(16) 

We leave the technical details to Appendix A. Important 
is that all terms are proportional to the stocked Casimir- 
Polder energy, that we will regularize as has been proposed 
in Eq. ©. The final result is, 

mvi = -P r a d = t]-(e - l)e inm Xmn (17) 
a 

with (see Appendix) r] = (J - h + C/3 - A/'i + D/3 - 
E/2)/l92ir 2 = 0.007909. 

2.1 Magneto-electric sphere 

We can insert the choice for a ME medium [22 , 23 : Xnm = 
gEu(E n B m - E m B n ) to find that 

mv = 2?7 x -(e - 1)3emE x B (18) 
a 

We recall that this relation applies after having turned on 
the fields E, B adiabatically. For a piece of FeGa03 (mass 
density 4.5 g/cm 3 , g^yiEB rss 10~ 4 ) of size a = 1/xm we 
find the unmeasurable speed v = 10 -20 m/s, some 12 
orders of magnitude smaller than the value predicted by 
Eq. UJ. 

2.2 Moving sphere 

It is well known that a sphere with a dielectric constant 
e, moving with a speed v much smaller than the speed of 
light, possesses a bi-anisotropic tensor Xij = (1 — £)cijk v k/co 
[4]. In this work we systematically neglect dispersion of 
the dielectric constant. We note however that the familiar 
dispersion for the dielectric constant beyond the plasma 
frequency, e = 1 — u) 2 P /uj 2 [27] will not be able to render 
the first Born approximation p5[) finite. This term would 
still diverge like J dkk. This suggests that the neglect of 
dispersion is not at the origin of the UV catastrophe, and 
that the problem is more fundamental. 

Equation (|17p thus applies and we can write, 

P rad = -2?7X A( e _i)2 v (19 ) 
ac 

The moving sphere thus drags along with him a radiative 
zero-point momentum with opposite sign. This is a rather 
revolutionary prediction, since it implies that the mass of 
the object is reduced by its finite size and its polarizabil- 
ity! We emphasize that this result follows from a dimen- 
sional regularization that is still arguably controversial. It 
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is difficult to say whether dispersion would eliminate the 
divergence. A scaling argument suggests that a dispersion 
as least as fast as e — 1 ~ 1/w 4 is required to give a finite 
result for the second order in the Born expansion. 

Note that in this particular case it is possible to book- 
keep the diverging momentum of the zero-point fluctua- 
tions into the bulk mass of sphere, i.e; the one measured 
for large a 3>. For a dielectric sphere of size 1 /ini the 



change is mass is of order 10 " 



and thus completely 



negligible. One could speculate about this effect on atomic 
scale (a = ao, a polarizability density of order so that 
e — 1 ~ 1 ). This would yield a mass reduction of order 
10 _4 m e (or 50 eV), still small enough to be unobserved. It 
is even more speculative to apply Eq. (fl~9|) on the electron 
scale (10 _15 m). Here the mass associated with zero-point 
momentum becomes of the same order of magnitude as the 
rest mass itself, that is of order m e = 0.5 MeV. Inspired 
by the original argument by Casimir to explain electron 
stability [15) one could even propose that the momentum 
of the electron is purely due to zero-point motion. If we 
propose P ra d = a(h/roCo)v := m e v we find that a should 
be equal to the fine structure constant e 2 /hco- Unfortu- 
nately, for the dielectric sphere we find the opposite sign, 
just like the Casimir force was also seen to be repulsive 
[5] , but it is fascinating that this argument gives the right 
order of magnitude for a. We remark that the value found 
above for r\ is close to the fine-structure constant. 



a(0) 



1 

a 



A, 



A similar "satisfactory regularization" procedure is not 
possible for the transverse Green's tensor. Dimensional 
regularization introduces r = 1/ At and the resonant fre- 
quency ui„ = c^a^r)- 1 / 2 . 

The regularization for the isotropic dipole produces 
an "acceptable" result for its scattering amplitude. In Ap- 
pendix B it is established that the inclusion of bi-anisotropic 
effects does not lead to new singularities. We can then use 
Eq. (fT4|) to calculate the vacuum expectation value for 
E x B. We find 



rad 



2?i T t 
— Im ^ v 

7T oalt 



(22) 



It can easily be checked that the frequency integral of 
Imto/iv 2 converges and equals — {-k /2)a{Q)uo/ c 2 . Hence 
we arrive at the final result, 



Pr 



ad 



a(0) hivo 



a 



(23) 



2.3 Moving dipole 

In the following we consider a moving electric dipole and 
calculate semi-classically the radiative momentum asso- 
ciated with the zero-point fluctuations, modified by the 
presence of the dipole. The polarization is assumed to be 
point like, and if the dipole is moving with speed v the 
light-matter interaction becomes 

V = (1 - e)^U\0) (0| - i-<t> P • x T U\0)(0\ 
eg c 

+z-(7|O)(O|x0 P (20) 

CO 

with U interpreted as a small physical volume associated 
with the dipole , a = (s — 1) its polarizability density, 
and the bi-anisotropic tensor \ = (1 — e)(e • v/c ). The 
advantage of this interaction is that the full Born series 
can be summed, although momentum integrals have to be 
regularized [28] . For v = 0, the t- matrix is found from 

i _ -A^r^/cl 

tQ ~ ~(cm; 2 )-i + G (r = 0) ~ u 2 - uj 2 - § iKfa/c~o 

(21) 

The second familiar formula is obtained when a mo- 
mentum regularization is adopted for the diverging k- 
integral of the Green's tensor [28]. The divergence of the 
longitudinal part Al/lo 2 > can be absorbed into the 
polarizability by defining 



We conclude that the moving dipole drags a momentum 
associated with zero point opposite to its kinematic mo- 
mentum. This can be interpreted as a reduction of the 
kinematic mass. The minus sign was also found earlier for 
the moving sphere. 

It is surprising to see that the ratio of real to bare po- 
larizability density comes in. For a small dielectric sphere 
difference between a and a(0) can be attributed to de- 
polarization induced by surface charges, but if we want 
to apply this the model to an atom, a is usually sup- 
posed to be unmeasurable. We can now imagine two sce- 
narios. If this depolarization is negligible, typically true 
when a ~ ct(0) rs U , the front factor in P ra( j equals one, 
and the mass would be reduced by an amount Hujo/cq. 
For a typical resonant transition at a few eV this would 
modify the Hydrogen mass by roughly one part in 10 9 . 
This is roughly the same value estimated in the previ- 
ous section on the basis of a dielectric atom. If however 
a 3> a(0) ~ U, the finite polarizability density is fully 
governed by surface depolarization, described here by the 
regularization scalar Al, then P ra d = 0. Unfortunately, 
the present semi-classical approach is not able to predict 
the value of a(0)/a. A quantum theory is needed. 

We expect in general that zero-point motion does not 
generate energy flow, not even in bi-anisotropic. This means 
that the quantum expectation value of the Poynting vec- 
tor should vanish. In a bi-anisotropic the latter is not nec- 
essarily proportional to the momentum. For the moving 
dipole it can be checked explicitly that the quantum ex- 
pectation value of the Poynting vector S = coE x H/4-7T, 
indeed vanishes. 
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3 Magneto-chiral object 

The optical properties of a homogeneous magneto-chiral 
(MC) material can be characterized by a contribution to 
the index of refraction that is independent on polarization, 
and linear in the magnetic field 29J Symmetry arguments 
impose that the sign of this contribution is different for 
opposite enantiomeres, and opposite for counter propagat- 
ing beams. Thus typically An ~ gh ■ Bo with g a material 
pseudo scalar related to microscopic chirality. This phe- 
nomenon can seen as a collective effect of rotatory power, 
with optical bi-anisotropy Xij = 9$ij j an< i the Faraday ef- 
fect that contributes iVeijkBk to the dielectric constant, 
with V the Verdet constant. 

If we accept this macroscopic description of optical 
MC, the "radiative" momentum of a MC sphere created 
by zero-point motion can be calculated in just the same 
way as was done earlier. The second order Born approx- 
imation generates optical MC by means of products of 
chiral and Faraday-type terms, that have to be regular- 
ized when we integrate over all frequencies of the vacuum. 
In a attempt to be more realistic one could accept that the 
Verdet constant behaves like V(ui) = Vouj 2 up to relatively 
large frequencies. On the other hand, rotatory power is ex- 
pected to have only little frequency dispersion. Without 
specifying details - the calculation is similar to the one 
for a magneto-electric sphere - we give here the outcome 
for the momentum obtained by a MC sphere, after having 
turned on the magnetic field adiabatically 



-0.005098 — ^B 



(24) 



If the assumptions and the regularization proposed 
above are correct, this formula would imply that enan- 
tiomeres of opposite chirality can be separated by turning 
on a magnetic field. 

A second approach consists of accepting the unavoid- 
able heterogeneous structure of space that underlies spa- 
tial chirality. One can propose a simple optical model to 
describe a chiral "molecule" in terms of a chiral distribu- 
tion of N > 4 classical dipoles. If the dipoles are subject to 
the Zeeman effect, this molecule exhibits MC properties in 
the optical scattering, that are particularly revealed when 
we average over orientations to restore spherical symmetry 

The scattering amplitude of the MC molecule - lin- 
earized in the external magnetic field - was obtained in 
Ref. 30J . It can be inserted into expression (|14p to find 
the radiative momentum. The end result can be expressed 
in terms of a trace of a complex 3N x 3iV matrix involv- 
ing two Levi-Civita tensor densities, and the scattering 
amplitude to of the dipoles found earlier. A straightfor- 
ward analyses leads us - quite surprisingly - to exactly 
the same expression as for the "diffuse supercurrent" that 
was considered by us in Ref. [3U]. The possibility of such 
a current, directed along the magnetic field and not in- 
volving the familiar gradient of energy density (familiar 
from Fick's law), was investigated for random media with 
chiral scatterers, but with negative result. In the present 



context we thus conclude that the chiral object does not 
carry any vacuum momentum, not at any frequency, when 
the magnetic field is turned on. 

One can try to analyse this conclusion. In the micro- 
scopic picture, Poynting vector and radiative momentum 
are proportional at any point, since locally /i = 1 and 
X = 0. Since we do not expect any macroscopic energy 
current - quantified by the average Poynting vector over 
some large volume - to occur in vacuum (yet this state- 
ment is hard to prove in heterogeneous, complex media), 
we might anticipate that also the macroscopic radiative 
momentum must vanish, 



Pr 



;ad 



47TCn 



|E x B|0) 



d 3 r(0|S|0> = 



(25) 



This second equality does not hold in the macroscopic 
description of MC, and the two deviate at any point. In 
the microscopic picture, dispersion and spatial structure 
have been taken into account much more realistically than 
in the macroscopic constitutive description. This example 
thus shows that one has to be careful in applying macro- 
scopic Maxwell equations to fundamental issues whose ori- 
gin is truly microscopic. The macroscopic, regularized out- 
come (IM1) is thus probably wrong. 



4 Conclusions 

The purpose of this work was to come to concrete ex- 
pressions for the momentum of zero-point motion in com- 
plex media. This constitutes a new and unique occasion 
to "test" regularization methods for vacuum properties 
in experiments, since momentum is directly observable, 
much more than energy. Three cases have been discussed 
for which the momentum of zero-point motion does not 
seem to vanish "trivially" . For a non-absorbing sphere 
subject to both an external electric and magnetic field 
we find a radiative momentum inversely proportional to 
its radius. To this end regularization techniques had to be 
adopted to render the outcome finite. A confrontation of 
this prediction to future experiments may thus shed new 
light on the validity of regularization methods in general. 
The same procedure leads to a radiative momentum of 
zero-point fluctuations of a moving sphere. This effect in 
principle lowers the kinetic mass of the sphere. The same 
conclusion is reached for a moving dipole, thus reassur- 
ing that - at least in this case - the prediction is not an 
artifact of the macroscopic model. At last, regularization 
techniques have been applied to a sphere exhibiting both 
rotatory power and the Faraday effect. Here it is possi- 
ble to come up with a more microscopic description, using 
Faraday-active dipoles in a chiral geometry. In this case 
the zero-point momentum is rigorously equal to zero, al- 
though the macroscopic models yields a finite result. 

In the futur we hope to develop fully quantum-mechanical 
descriptions of magneto-electric objects and moving dipoles. 
The calculations have also been done for idealized media, 
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free from dispersion and absorption. Clearly, this has to 
be improved in the future. The Lorentz-invariance of zero- 
point motion is also an important aspect that must be 
given attention. 

The author is indebted to Geert Rikken for many enlightening 
discussions. This work was supported by the European Ad- 
vanced Concept Team, call ACT-RPT-ARIADNA-04-1201. 
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A Calculation of Iyi 

In this Appendix we calculate the iyj defined by Eq. (JT5 
For simplicity we put cq = 1. Since the constitutive pa- 
rameters are assumed not to vary with frequency, the 
frequency integral of this object can be performed using 
Cauchy contour integration. Since Go = (1 — kk/w 2 )(a; 2 — 
k + ie) has a longitudinal part and a part proportional 
to the identity, we can essentially discriminate three dif- 
ferent contributions to Eq. fTB")) , with either 0, 1, or 2 
longitudinal propagators in the above expression. 

In the absence of any longitudinal part, the frequency 
integral of any contribution proportional to (s — 1)\ con- 
tained in (fT6"| is of the form 



rfww 4 



1 



1 



iri k + 2k' 

/(l ( w 2 _ k 2 + ie f uj 2 - k' 2 + fe ~ ~T (k + k') 2 

(26) 

We can write V(w, k',k) = V(k, k')#kk' with 



#kk' = / d 3 xexp[i(k — k') • x] 
Jb 



(27) 



where in our case the integral is to be carried out over a 
sphere. This transforms Eq. (fT6|) into, 



I 



(0) 



1 



v l 4 



(e- l)(e jmn Xlm +Xjm^lmn) / d x / d y 



/ (2^3 / (4^3 ex p[* k • ( x - y)] cx pH k ' • ( x - y)] 



k + 2k' 



The four integrals restore complete spherical symmetry, so 
that the tensor ki(k n + k' n ) must lead to a factor propor- 
tional to 5i ra /3. The factor of proportionality then imme- 
diately follows by contraction. If we re-scale fc = |x — y\p 
and k' = |x — y\q we arrive at the following expression, 

T ijl = ""^"X^ ~ 1 )( io ~ Jl)( e im»Xim + Xjmei mi )(28) 

where we have introduced the volume integral 

K{B) ee f d 3 x f d 3 yr— ^ 
Jb Jb l*-y| 7 

and the two scalars 



Io= l°° dp [°° dqpV 7~~T~~T2 Jo (P)J0 (<?) (29) 
Jo Jo (P + q) 

h = J a dp dqp 3 q 3 ji(p)jt(q) (30) 

A factor Q~ e ^P+i) can be added to ensure convergence. 

We can repeat this calculation in the presence of one 
longitudinal propagator (— kk/o> 2 )(ijj 2 — k 2 + ie)^ 1 . The 
frequency integral now becomes, 



du> uj 2 



1 



iri 1 



{uj 2 - k 2 + ief uj 2 -k' 2 + ie 4 k(k + k') 2 

(31) 
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We give the end-result of a long calculation that involves = 8.246 • • • 
angular averaging of four-rank tensors, that generate Bessel r™/ 2 

functions of higher order, and that make the calculation E 2 = -4 / cZ0cos 2 0(3 + 1.5 sin 20 sin 40 

of scalars more involved. We find, 

+ 3sin0sin30 + cos0cos30) = —13.744 ■ 



L iji I67r 4 E 3 = 6 c?0cos0(6cos0 — 2cos20cos30 







C — A\ 

(Di H f^ - J Xmn(8ije n lm + Siitnjm) - sin 2 20cos50) = 24.74- 

C-A\ 

D3 H — — I [emliXjm + ^mjiXlra) 



(32) 

B Radiative momentum of moving dipole 



in terms* of D 'Andl2Di + ?>D 3 E E\ + 2 + E3 The scattering matrix of the moving dipole can be calcu- 

1 lated from the Born series (c = 1) 

f°° f°° p 4 q 3 
A= / dp dq- — ■ — rnji{p)ji{q) 

Jo Jo (p + qr r d 3 k" 

r°° r 00 rPn 2 1W = V kk < + / „ V kk » ■ Go(k") ■ V k » k > H — 

C = / ^ / j (p)j ( g ) J ^ 



(p + q) 2 



which for the bi-anisotropic point dipole can be fully summed 



f°° f°° p 3 q 4 
D = J ^ J ^ (p + j)^ ' ™^ UP t0 ' 

E = 1 / d 3 /■ d 3 1 (P^O! [i( } . f] T kk , = i + — ^- [-z0 k • X T + z X k ,] (36) 

(33) 

Here, io is the i-matrix of the isotropic dipole, given in 

™ ,. , , . , , j- 1 the text. We can develop Eq. (fl4l) to 

Ine final contribution involves two longitudinal prop- ' — ' 

agators. Since the angular average over a tensor of order 6 

appears it is convenient to perform the contractions first. fa r ^3^. 

We find that, Q ' PradM = --ImTr J ^—3 (k ■ q - kq) ■ G (k) ■ 
(2) _ (2) _ K(B)E [(e . k) • (e • v) + (e • v) • (e • k)] ■ G (k) 

ij'j J»i ~~ 327T 4 JtimnXnm V J V 

We have (e • k) • (e • v) + (e • v) • (e • k) = kv + vk - 2k • v. 
We conclude that all contributions generate the stocked It is convenient to use the identity, valid to order v, 
Casimir-Polder energy K{B) of Eq. @ as the only diverg- 
ing element. If we adopt dimensional regularization for ^ s r , , . „, i <-< n \ t~> n . ^ ^ n ^ \ 
thl object, G (k).[-kv - vk + 2k • v].G (k) = G„(k+-v)-G (k--v) 

7T 2 

^ ~ * 12a ^" n S s us to 

To calculate the double momentum integrals above we % r rf 3 k 

can move the p- integral to the imaginary axis and formu- 1 ' Prad(^) = ~I m Tr J 7^33 (l ' v — Q v ) ' G*o(k) 
late the resulting integral in the full complex plane using 

polar coordinates (r, 0). We find Since j d 3 k ( 2vr )-3 Go (k) = -\juu? -l/t this reduces to 

Io = -12 / ' d0 cos 50 cos sin 4 = 0.589- •• q . P rad US) = (q . v) x — Im (37) 

Jo 

f7r/2 



Ji = — 6 / <i0 sin 2 cos 0(3 sin 20 sin 50 + 2 cos 30) 
Jo 

= 4.123- ■■ 

A = 4 dxj) sin 2 cos cos 20(2 cos 30 + 3 sin 20 sin 50) 
Jo 

= 1.374- • 

C = 24 / d0sin 4 0cos 0cos 50 cos 20 = — 1.767 • • • 
Jo 

,tt/2 

i?i = — 24 / <i0 sin 2 cos 0(1.5 sin 20 sin 50 + cos 30) 
Jo 



